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ABSTRACT 


The  physical  model  of  a  multi-constant  elastic 
material  [1]\  with  consideration  of ' couple-stresses, 
is  studied  critically  to  determine  whether  it  more 
reasonably  represents  a  3-constant  isotropic  material 
or  a  6-constant  centro- symmetric  cubic  crystalline 
material.  A  description  of  a  critical  shear  experiment 
on  the  model  is  presented  along  with  experimental  re¬ 
sults.  A  parameter  a  which  is  definitively  related 
to  the  question  is  defined  and  ifcs  value  determined  on 
the  basis  of  statical  deformation  experiments.  The 
conclusion  is  that  the  model  behaves  as  a  cubic  crystal 
showing  couple- stress  effects. 
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Numbers  in  brackets  refer  to  references  at  end  of 
report . 


INTRODUCTION 


In  a  previous  paper  [l],  a  physical  model  of  a  so-called 
3-constant  isotropic  elastic  material  was  used  to  illustrate 
the  phenomenon  of  couple- stresses.  The  model  and  its  analysis 
were  predicated  on  the  assumption  that  couple-stresses  exist 
and  that  the  model  represents  an  isotropic  material,  at  least 
to  some  approximation.  An  obvious  question  provoked  by  the 
structured  form  of  the  model  is  whether  it  does  not  more 
accurately  represent  a  cubic  elastic  material.  As  previously 
pointed  out  [l].  Professor  R.  D.  Mindlin  strongly  urged  con¬ 
sideration  of  such  a  possibility.  Consequently  the  present 
paper  is  intended  to  provide  an  analysis  of  the  model  from 
such  a  viewpoint.  At  first  it  is  considered  worthwhile  to 
briefly  review  the  classical  theory  of  the  constitutive  equa¬ 
tions  of  anisotropic  materials  and  its  extension  to  include 
considerations  of  couple-stresses. 


CONSTITUTIVE  EQUATIONS  OF  GENERAL  ANISOTROPIC  SOLIDS 

The  classical  theory  of  the  deformation  of  crystals, 
without  consideration  of  couple- stresses,  has  been  extensively 
studied  in  the  past  [2,  3,  4],  The  relationship  of  the  con¬ 
cept  of  Isotropy  to  that  of  a  crystalline  structure  is  a 
peculiar  one.  The  concept  of  an  isotropic  body  is  not  con¬ 
tained  in  that  of  the  crystalline  group  and  can  only  be 
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defined  as  a  body  resulting  from  a  random  orientation  of 
crystals  which  produces  a  statistical  equivalent  to  a  solid 
continuum  specified  by  two  independent  elastic  constants.  It 
is  clear  that  if  one  contemplates  the  perfectly  crystalline 
structure  he  diverges  from  the  concept  of  isotropy.  The 
group  of  crystals. extends  from  the  simplest,  which  is  the 
centro- symmetric  cubic,  to  the  most  complex,  which  is  char¬ 
acterized  by  21  elastic  constants  [4].  The  defining  device 
for  crystalline  structure  is  found  in  the  mathematical  trans¬ 
formations  of  rotation,  reflection,  and  inversion  [4].  In 
the  classical  theory,  the  constitutive  equations  relating  the 
strain  tensor  and  the  stress  tensor  contain  8l  constants 
which,  because  of  symmetries  and  energy  considerations,  re¬ 
duce  to  21  constants  in  the  case  of  the  most  complex  crystal 
[2].  For  the  simplest  crystal,  the  centro- symmetric  cubical 
one,  the  independent  constants  reduce  to  3  in  number.  As  is 
well-known,  there  are  two  independent  constants  in  the  case 
of  isotropic  materials. 

As  long  as  couple- stresses  are  considered  to  be  non¬ 
existent,  the  old  classical  theory  of  crystals  is  probably 
satisfactory  for  most  purposes.  However,  once  couple- stresses 
are  postulated,  as  t^e  present  authors  consider  that  they 

”  “ 

should  be,  the  situation  becomes  more  complicated.  Kroner 
assumes  the  burden  of  this  complexity  in  his  attempt  to  ex¬ 
plain  deformation  phenomena.  He  uses  the  idea  of  the  motion 
of  dislocations  associated  with  the  assumed  existence  of 
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couple- stresses  to  explain  the  phenomenon  of  plasticity. 
Furthermore,  he  also  speculates  that  the  deformation  of  the 
elastic  framework  of  perfect  crystals  may  introduce  the  need 
for  couple- stresses  in  order  to  maintain  equilibrium  [5]* 

The  view  of  the  present  authors,  however,  is  that  a  student 
of  continuum  theory  is  not  constrained  to  espouse  any  partic¬ 
ular  theory  of  the  mechanical  structure  of  matter  but  may  at 
once  resort  to  the  generalization  of  the  concept  of  continua 
by  postulating  the  existence  of  couple-BtresBes.  From  this 
point  of  view  the  onus  is  on  anyone  who  denies  the  existence 
of  couple- stresses  to  demonstrate  that  they  do  not  exist  in 
any  particular  deformed  body. 

In  any  event  the  present  paper  is  concerned  with  the 
proper  definition  of  a  cubic  elastic  material,  in  the  presence 
of  couple- stresses  and  its  analysis.  Furthermore,  having 
determined  the  explicit  form  of  the  constitutive  equations 
for  such  material  the  deformation  of  the  physical  model  [1] 
will  be  studied  to  determine  whether  it  approximates  a  cubic 
crystal  in  itB  action. 


DERIVATION  OF  THE  CONSTITUTIVE  EQUATIONS 
FOR  A  6-CONSTANT  CUBIC  CRYSTAL 

For  the  purpose  of  clarity,  a  few  words  about  the  defini¬ 
tion  of  stress  and  strain  will  now  be  given.  Because  the 
existence  of  couple- stresses  is  postulated  in  the  theory,  it 
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will  be  necessary  to  use  symbols  for  the  symmetric  part  and 
the  anti- symmetric  part  of  the  stress  tensor.  Furthermore, 
because  of  the  peculiar  nature  of  couple- stresses,  the  couple- 
stress  tensor  will  be  divided  into  the  deviatoric  and  the 
spherical  parts  [6]. 

The  usual  total  resultant  shearing  and  normal  stresses 
will  be 'designated  by  t . .  and  written  in  terms  of  the  sym- 
metric  arid  anti- symmetric  parts  as  follows: 


(1) 


where  is  the  symmetric  part  and  is  the  anti¬ 

symmetric  part. 

Furthermore,  the  couple- stresses  will  be  designated  by 
qij  and  written  in  terms  of  the  deviatoric  part  and  the 
spherical  part  as  follows: 


=  q 


i,1  +  3  (qkk)6ij 


(2) 


where  q.  is  the  deviatoric  part  and  the  spherical  part  is 
represented  by  the  terms  in  .  The  repeated  index  means 
summation  over  the  3- dimensional  indices  1,  2,  3  >  as  usual. 


and  6 


ij 


is  the  Kron^-cker  delta. 


Now  ^  TjA  in  general. 
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It  may  be  observed  that 


+ 


and 


U 


(3) 


The  lineal  and  shearing  strains,  having  the  same  meaning 
as  in  classical  elasticity  [7]*  will  be  designated  by  e 
and  are  written  in  terms  of  the  displacements  u^  as  follows: 


+ 


u 


w 


where  the  comma  means  differentiation.  The  deformation  corre¬ 
sponding  to  couple- stresses  will  be  designated  by  and 

written  as  gradients  of  rotation  as  follows: 


Kij  eimj  “mi,! 
where 


(5) 


0  ,  when  any  two  of 

+1  ,  when  i,m,j  is 
numbers  1,2,3 


=  -1  ,  when  i,m,j  is 
numbers  1,2,3 


the  indices  are  equal 
an  even  permutation  of  the 

an  odd  permutation  of  the 


“W  -  f(\j  -  ui,m> 

or  as  some  investigators  do  we  may  write 
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Um 


(6) 


) 


KU  =  ®J,i 

where  in  rectangular  Cartesian  coordinates  the  components  of 
rotation  are  as  follows: 


1 — 1 

au^ 

2 

1 

au, 

(~i 

3x3 

CVJ 

1 

au0 

2 

aur 

_ C. 

ax- 

dx- 


Su- 


bxr 


(7) 


For  a  centro- symmetric  general  anisotropic  material  the 
energy  function  W  may  be  written  for  an  initially  unstrained 
elastic  body  as  follows: 


-  V 


_1_ 

2 


c  pk-ft  i  1  rpf  J  V- 

Ski  €i j  G  +  T  TJJ m  K1  j  K 


(8) 
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where  and  T ^  are  elastic  constants  corresponding  to 

the  force-stresses  and  the  couple-stresses  respectively  L 6 ] . 


From  the  energy  function  expressions  for  the  symmetric 
part  of  the  stress  components  and  the  deviatoric  part  of 

j 

the  couple-stresses  may  be  written  as  follows: 


aw 


'ij 


and 


U 

aw 


ax 


M  . 
SiJ  kl 


4m 

Tij  Km 


(9) 
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In  a  rectangular  Cartesian  system  of  coordinates: 


(10) 


Also  the  usual  symmetry  conditions  can  be  proved  [3]: 


,ij  _ 


TU  =  T^m 


(11) 


The  4th  order  tensors  S~^  and  in  indices  which 

run  over  the  values  1,2,3  for  3- space  have  apparently  8l 
components.  However,  as  is  well-known  symmetry  and  thermo¬ 
dynamic  conditions  reduce  the  number  radically.  In  fact,  it 
is  readily  shown  that  the  8l  values  for  reduce  to  21  by 

the  symmetry  relations  [7]  and  the  further  symmetry  relations: 


(12) 


which  result  from  the  symmetry  relations: 


s  _  s 
Tij  “  TJi 


€ij  " 


(13) 
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The  further  reduction  of  the  tensor  to  represent 
the  case  of  the  centro- symmetric  cubic  crystal  is  given  in 
standard  texts  on  crystallography  [4].  Then  the  three  inde 
pendent  nonzero  components  turn  out  to  be: 


11 

S22 

_  o33 

11  “ 

b22 

“  S33 

23  _ 

s13 

-  S12 

’23  “ 

b13 

b12 

,11 

ell 

c22 

’22  “ 

S33 

"  S33 

(14) 


Using  contracted  notation  for  simplicity  of  writing  the 
explicit  form  of  that  portion  of  the  constitutive  equations 
which  derive  only  from  forces  and  not  from  couples,  may  be 
written: 


lll  " 

°11 

11 

T 

°12 

22 

°12 

33 

ts  = 
t22 

S12 

6 11 

+ 

S11 

€22  + 

S12 

e33 

tS  = 
t33 

S12 

€11 

+ 

S12 

€22  + 

S11 

€33 

ts  = 

t23 

S44 

e23 

t13  = 

S44 

€13 

T12  = 

S44 

£12 

(15) 


Now  to  derive  that  portion  of  the  constitutive  equations 
which  derive  only  from  consideration  of  couples  at  a  point, 
it  is  necessary  to  take  into  account  the  asymmetry  of  K. . 
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and  q!j\j  as  well  as  the  geometrical  fact  that: 

Kii  -  0  (16) 

which  implies  that 

lii  =  o  .  (17) 

[again  repeated  index  implies  summation]. 

It  is  then  clear  that  the  reduction  of  the  4th  order  ten¬ 
sor  Tj^  is  formally  somewhat  different  from  the  reduction  of 
Sj^  .  From  the  conditions  (l6)  and  (17)  on  and  re¬ 
spectively  the  following  nine  conditions  are  Imposed  on  : 

ki 

=  0  (18) 

where  i,j  take  on  the  values  1,2,3  and  the  repeated  index 
k  implies  summation. 

From  the  methodical  application  of  the  transformations 
of  rotation,  reflection,  and  inversion,  the  explicit  form  of 
that  portion  of  the  constitutive  equations  which  derives  from 
consideration  of  couples  at  a  point  may  be  written  as  follows: 
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d 

qll 

= 

3A  v- 

7  Kn 

_d 

q22 
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3A  «■ 

7  *22 

_d 

q33 
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3A  K 

7  K33 

d 

q23 

— 

+  BK^2 

_d 

q32 

= 

CK32  +  BK^ 

_d 

q13 

= 

ck13  +  bk31 

q31 

= 

CK31  +  BK1 3 

<4 

= 

CK12  +  BK21 

oil 

= 

CK21  +  BK12 

(19) 


where  A,  B,  C  are  independent  constants. 

If  the  material  were  isotropic,  the  following  conditions 
apply: 

S11  "  S12  “  <S44 

and  (20) 

^  =  B  f  C 
2 

For  a  centro- symmetric  cubic  crystal  the  six  independent 
elastic  constants  are  then  S^1  ,  S12  »  ,  A,  B,  and  C  . 
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As  a  first  step  beyond  the  theory  of  classical  crystal¬ 
lography  which  requires  only  the  constants  ,  S12  ,  , 

it  seems  appropriate  to  place  B  equal  to  zero.  In  this 
event,  the  number  of  nonzero  independent  constants  is  appar¬ 
ently  five,  which  are  the  three  corresponding  to  classical 
elasticity  and  the  two  introduced  because  of  consideration  of 
couples  at  a  point.  These  latter  two  constants,  A  and  C  , 
correspond  to  twisting  couples  and  flexure  couples  respec¬ 
tively. 


PHYSICAL  MODEL  OF  A  6-CONSTANT  CUBIC  ELASTIC  MATERIAL 


As  observed  in  the  introduction  to  the  present  paper, 
one  of  its  principal  objectives  is  to  examine  the  question 
of  whether  a  physical  model  used  previously  [l]  represents  a 
3-constant  isotropic  material  really  or  more  appropriately  a 
6-constant  cubic  elastic  material.  A  crucial  test  consists 
of  a  determination  of  the  value  of  a  parameter  a  given  as 
follows  [9l: 


S  44 

S11  '  S12 


(21) 


If  a  equals  unity,  the  first  condition  (20)  is  satisfied 
and  the  material  is  presumed  to  be  isotropic.  On  the  other 
hand  if  a  is  not  equal  to  unity,  with  due  consideration  of 
margins  of  experimental  error,  the  material  will  be  presumed 
to  be  cubic  crystalline. 
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As  it  turns  out,  the  tension  experiments  previously  de¬ 
vised  and  conducted  [1]  properly  determine  the  values  of  the 
constants  ,  S^g  whether  the  body  is  isotropic  or  cubic. 

Therefore  only  a  single  additional  experiment,  for  the  deter¬ 
mination  of  the  shear  modulus  ,  is  required.  The  dis¬ 

placement  functions  and  the  shear  experiment  will  now  be 
described.  For  convenient  reference  some  details  of  the 
physical  model  are  shown  in  Fig.  1. 

EXPERIMENT  TO  DETERMINE  ELASTIC  SHEAR  CONSTANT  S^ 

In  order  to  evaluate  the  shear  constant  S^  it  is 
sufficient  to  assume  that  a  rectangular  prism  is  loaded 
around  its  lateral  faces  with  uniformly  distributed  shearing 
stresses  of  constant  magnitude,  say  t  as  shown  in  Fig.  2. 

The  displacements  caused  by  such  a  loading  may  be  obtained 
theoretically  by  solution  of  the  corresponding  plane  strain 
problem  of  an  infinitely  long  rectangular  prism.  It  is  as¬ 
sumed  satisfactory  for  experimental  purposes,  to  an  acceptable 
degree  of  approximation,  to  represent  such  a  theoretical 
situation  by  a  slice  of  the  prism. 

In  order  to  obtain  theoretically  determined  displacement 
functions  it  is  necessary  to  solve  the  six  equations  of  equi¬ 
librium  [1,  81,  three  on  force  and  three  on  moments.  Obviously 
for  the  present  assumed  shape  and  simple  loading,  these  equa¬ 
tions  reduce  to  the  following  two: 
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0 


Now,  as  usual,  an  integration  of  the  strain- displacement 
equations  will  provide  the  necessary  equations  for  the  deter¬ 
mination  of  the  shear  modulus  .  The  equations  are: 
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du. 

dx- 

du, 

< 

dx,. 


'2 


and  a  &  +  !^) 

2  dx2  “dx^ 


544 


(25) 


Consequently  one  obtains  for  the  displacement  functions  u.^ 
and  u2  the  followings 


and  u2 


(26) 


The  displacement  (u1  ,  u2)  can  be  measured  at  any  point 
(x^^  ,  x2)  with  the  model  loaded  by  constant  shearing  stress 
t  .  Consequently,  can  then  be  evaluated. 

A  model  consisting  of  metallic  cubes  connected  by  either 
flat  steel  strips  or  round  steel  rods  was  designed  and  fabri¬ 
cated.  It  is  of  the  same  type  as  the  one  previously  described 
in  the  literature  [l].  For  the  present  investigation,  two 
types  of  connector  pieces,  flat  and  round,  were  used.  The 
reason  for  making  model  slices  with  the  round  rods  is  that  it 
is  planned  to  build  up  later  large  prisms  which  are  essentially 
3-dimensional  models  for  study  purposes.  Such  models  will 
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enable  one  to  measure  the  final  elastic  constant,  the  twist 
modulus  A  ,  which  is  still  undetermined. 

As  one  may  recall  from  the  previous  paper  [l],  the  flat 
steel  connector  pieces  had  cross-  sectional  dimensions 
0.5U  by  0.016".  The  round  steel  rods  used  in  the  second 
model  are  0.095 ‘inches  in  diameter  and  2  inches  between 
cubes. 

In  order  to  perform  the  experiments  each  model  was 
mounted  horizontally  so  that  it  rested  on  several  steel  ball 
bearings  in  order  to  maintain  itself  in  the  horizontal  plane 
when  it  was  loaded  and  at  the  same  time  minimize  any  frictional 
forces  which  might  oppose  the  free  deformation. 

In  order  to  deform  the  model  as  required,  loads  were 
applied  through  special  fittings  such  that  they  were  uniformly 
distributed  through  the  thickness  and  along  the  sides.  The 
experimental  apparatus  and  model  are  shown  in  Fig.  3.  The 
elongation  and  contraction  of  the  two  diagonals  were  measured 
using  special  displacement  gages  [1].  Elongation  of  the  one 
diagonal  and  simultaneous  contraction  of  the  other  were  prac¬ 
tically  equal  in  magnitude  as  expected.  However,  on  account 
of  small  differences  the  two  measured  lengths  were  averaged  in 
order  to  calculate  the  change  in  angle  7  at  any  corner  of 
the  model.  Measured  values  are  shown  in  Fig.  4.  For  infini¬ 
tesimal  deformation  theory,  the  change  in  angle  is  given  by 
the  expression: 
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(27) 


3u.  3iu 

(_1  +  — — ) 

dx2  dx^ 

Hence  knowing  the  value  of  the  change  of  angle,  the  value 
of  the  constant  S^  can  be  readily  calculated  from  Eq.  (26). 

An  additional  experiment  was  conducted  on  the  model  with 
round  rod  connectors  in  order  to  determine  the  values  of  S-^ 
and  S12  . 

All  of  the  elastic  constants  so  far  determined  by  means 
of  experiments  are  shown  in  Table  I. 


DISCUSSION  AND  CONCLUSIONS 

Since  the  pertinent  .elastic  constants  have  been  determined, 
the  decisive  constant  a  needed  to  determine  whether  the  model 
may  be  considered  isotropic  or  cubic  crystalline  can  be  calcu- 

-TJ 

lated.  Its  value  for  the  first  model  is  19  x  10  and  for 

—  ^ 

the  second  model  15*8  x  10~“>  . 

It  is  clear  that  since  the  values  of  a  are  so  radically 
different  from  unity,  the  models  cannot  be  considered  as  repre¬ 
senting  an  isotropic  material.  So  far  as  the  evidence  goes 
they  may  be  considered  as  models  of  the  somewhat  more  complex 
cubic  crystalline  matter. 

It  should  be  recalled  that  the  value  of  the  constant  S12 
for  the  present  type  of  model  is  approximately  zero. 
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The  last  constant  associated  with  a  cubic  crystalline 
material  exhibiting  the  phenomenon  of  couple- stresses  is  the 
twist  constant  which  has  been  designated  A  in  the  present 
paper.  For  the  purpose  of  determining  the  value  of  this  con¬ 
stant  it  will  be  necessary  to  fabricate  a  3-dimensional  model 
built  up  out  of  t*he  cubelike  substructures.  Such  a  model,  of 
sufficiently  large  proportions,  is  now  being  constructed  and 
probably  will  be  subjected  to  an  experiment  in  the  near  future. 

In  addition  to  experiments  for  the  determination  of  the 
twist  modulus  A  ,  vibration  and  elastic  wave  transmission 
experiments  are  planned. 
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Table  I.  Elastic  Constants 


Model 

Tensile 

Modulus 

Shear 

Modulus 

Flexure 

Modulus 

S^h  lbs. /in. 

S^h  lbs. /in. 

Ch  lbs. /in. 

With  steel  strips 

as  interconnection 

.182  x  106 

3572 

.295  X  105 

With  steel  rods 

as  interconnection 

.132  X  106 

2080 

.777  x  105 

10  x  10  CUBE  MODEL 


Some  Details  of  Physical  Model 


Fig.  2 .  Section  of  Model  Showing 
Type  of  Loading 


dBBB 


(•sqx)  sspis  aqq  uo  p^oq  tbioj; 


